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Abstract. The recent development of numerical schemes for Relativistic 
MHD (RMHD) allows us to model the acceleration and outflow properties of 
winds from compact sources. Theoretical models suggest that acceleration and 
collimation of the flow are extremely inefficient when the speed is close to c, in 
contrast with many observations. Numerical results for an axisymmetric rotator, 
both in the case of monopolar and dipolar magnetic field will be presented, sug- 
gesting that in ideal RMHD acceleration is indeed inefficient. I will also mention 
numerical challenges and stability problems of present simulations. Finally I will 
discuss the interaction of a pulsar wind with the surrounding SNR. I will present 
emission maps based on numerical simulations of the flow inside the nebula and 
how they can be used in comparison with observations to derive informations 
about the properties of pulsar winds far away from the source. 



1. Introduction 

Magnetically dominated outflows from stars and accretion discs, are central to 
the evolution and properties of these objects. Because they can efficiently extract 
rotational energy, and convert it into wind kinetic energy, magnetic outflows 
are ubiquitous in powering a variety of astrophysical systems. Rotationally 
powered pulsars are among the best example of relativistic outflows produced 
by a fast compact rotator. The high magnetic field at the surface (B ~ 10 10 ~ 13 
G) and short rotation period (P ~ 0.001 — Is), induce an electric field at the 
surface itself which is strong enough to lift electrons and ions and accelerate 
them to relativistic speeds. These accelerated particles, as they move along 
curved magnetic field lines, produce the pulsed emission that is observed. In the 
magnetosphere of the neutron star, these particles produce pairs. If the charge 
density exceeds the Goldreich-Julian value then conditions can be modeled in 
term of ideal RMHD. 

The typical Lorentz factor of the plasma in the pulsar magnetosphere, 
within the light cylinder, is 7 ~ 100. However from observations of the in- 
teraction of the pulsar wind with the environment, and especially the surround- 
ing SNR, we know that far away fr om the source the wind is ul trarelativisti c 
(7 ~ 10 4-6 ) and weakly mag netized (|Kennel fc Cornonitilll984aH bl: lAronsll200^ . 



The basic idea is that the wind is magneto-centrifugally accelerated, and that 
the high Lorentz factor is achieved by almost completely converting magnetic 
energy into kinetic energy. 

In this sense, pulsars constitute the prototype of relativistic accelerators 
and the basic physics of magneto-centrifugally driven relativistic winds apply 
to a wide variety of astrophysical objects from accretion disks to GRBs and 
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microquasars. In all of these systems a plasma is assumed at the base with small 
outflow velocity (compared to the asymptotic velocity), embedded in a strong 
magnetic field. Magnetorotation then accelerates the plasma and converts the 
magnetic energy into relativistic motion. The study of winds from neutron stars 
can give an useful insight into the general problem of relativistic outflows. 

In the simp le ID radial model for relativistic winds from compact rotators 
hv lMichell (|l969T l. solutions can be parametrized in term of a = S1 2 <J> 2 /M, where 
Q is the rotation rate, <3? the magnetic flux and M the mass flux, a is the 
maximum Lorentz factor the wind can achieve if all magnetic energy is co nverted 
into k inetic energy, and at the base of a pulsar it is ~ 10 6 . The results bv lMichell 
(1969) show that, due to relativistic effects, the asymptotic Lorentz factor tends 
to er 1 / 3 , and conversio n of magneti c to kinetic energy is minimal. This is known 
as the a — 7 problem (|Aronsll2004h . 

The low efficiency of the Michel model is due to the assumption of strictly 
radial flow. If this assumption is relaxed and the flow is allowed to diverge 
more than radi ally from the Light C ylinder up to the SNR, then acceleration is 
more efficient ((Begehnan Sz Lilll994h . It is known from the theory of classical 
winds that the presence of a magnetic field can cause collimation of the flow 
along the rotation axis, and divergence along the equator, so in principle one 
might expect magnetic collimation to be also responsible for the acceleration of 
the wind. However when the flow speed reaches values close to c the Coulomb 
term in the electromagnetic force cannot be neglected (as in non relativistic 
MHD). The strength of the induced electric field is comparable with the magnetic 
field and the Coulomb term balances the magnetic hoop stress. Solutions of 
the rela tivistic Grad-Shafra nov equation for a monopolar field in the trans-fast 
regime ( Beskin et al. I flQ98T > have shown that acceleration is too slow and the 
asymptotic results of the Michel model apply also in 2D axisymmetric geometry, 

meaning that the flow remain radial. 

All simulations were performed using the code developed bv lDel Zanna et alJ 



(2003). It is a high-order finite difference UCT scheme, that solves the equation 
of relativistic MHD in general curvilinear coordinates. The scheme is based on 
an approximate HLL solver, that avoids wave decomposition. 



Wind acceleration 



The first axisymmetric n umerical simulatio ns of the acceleration of a relativistic 
wind were presented by iBogovalovl ((1999), where the case of a zero pressure 



plasma in a monopolar magnetic field, with uniform injection conditions at the 
base of the star was investigated. Results show that when 7 exceeds 5 collimation 

and acceleration of the wind become marginal. 

Our recent investigation ( Bucciantini et af1l2006h has focused on the study 



of the acceleration of relativistic wind in the case the pressure is small compared 
with the rest mass energy density but not negligible. It has been proposed that 
soon after the birth of a proto neutron star (NS), neutrino heat ing might rise the 
soun d speed in the neutron star atmosphere to value ~ 0.1c (|Thompson et all 
2004). The heating might be able to drive a wind against the gravitational at- 
traction of the NS. Such wind might be accelerated to high 7 by the magnetic 
field. If acceleration is efficient, the wind can reaches 7 ~ 100 — 300, enough to 
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Figure 1. Wind structure for a monopolar magnetic field, a — 200, axes are 
in units of r/vs- Left: contours represent magnetic surfaces, colors represent 
{B^/Brl, the dashed line is the Alfvenic surface, coincident with the light 
cylinder, the dotted line is the fast-magnetoson ic surface. Right: Lorentz 
factor in the wind. From lBucciantini et alJ ((2006) 



trigger a GRB. We have investigated the axisymmetric wind from rotating neu- 
tron star in Schwartzschield metric. Pressure, density and the radial component 
of the magnetic field are fixed at the NS surface (r^s)- A 2D spherical r — 
grid was used, with 100 uniform cells in 9, and a logarithmic scale with 200 cells 
per decade in r. 

Solutions in the 2D case can be parametrized by using the value a averaged 
over all stream lines. Results in the case of a monopolar magnetic field show 
that for a > 1 the flow appears to be radial, the terminal 7 is <C a, pointing 
to inefficient acceleration. Moreover the energy and angular momentum fluxes 
rapidly approach the force-free value. Despite the fact that in the close mag- 
netosphere v <C c, asymptotic results agree with the flow structure found by 
iBogovalovl Jl999). Energy flux scales as sin (#), and 7 oc sin(#), in agreement 
with the exact monopole solution. Results suggest that a rapidly rotating mag- 
netar will not be able to simultaneously produce a collimated and relativistic 
outflow. Relativistic velocities can be reached only for £r>l, and in this case 
most of the energy would be concentrated in the equatorial plane instead that 
in a jet. 

We computed solutions also in the more realistic cases of a dipolar magnetic 
field at the surface of the NS. Due to stability issue in numerical RMHD it was 
not possible to study cases with a > 20. Numerical resistivity on current sheets 
at the edge of the closed zone, causes the plasma to rotate faster that corotation, 
while at the equator reconnection takes place and plasmoids are formed. To 
obtain steady solutions one is forced to suppress resistivity, thus enhancing the 
instability of the code. Moreover, in the case of a dipole the magnetic and flow 
surface are not aligned with the grid, and in this case inversion from conservative 
to primitive variables is not stable for B 2 /pc 2 > 100. 
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Figure 2. Wind structure for a dipolar magnetic field, a — 20, axes are 
in units of r^g. Left: contours represent magnetic surfaces, colors represent 
\B<p/B r \, the dashed line is the Alfvenic surface, coincident with the light 
cylinder, the dotted line is the slow-magnetosonic surface. Right: same as 
the left panel, the dotted line is now the fast-m agnetosonic surface . Note the 
radial structure in the asymptotic region. From .Bu cciantini et alJ l)2006j) . 



Despite the presence of a closed zone, and the fact that within the light 
cylinder the field structure is close to a dipole, we found that solutions can 
again be parmaterized in term of a, now denned as an average done only on 
open streamlines. In this case we find that energy and angular momentum 
losses are the same of a monopolar case wit h the same a. Again as soon as 
a > 1 solutions approach the force- free limit ( Contopoulos et al1ll999h . outside 



the Light Cylinder, and for a > 10 the asymptotic structure resembles the 
exact monopole, with energy flux peaking at the equator and a highe r equatorial 
Lorentz factor. Similar results are found also bv iKomissarovl (|2006f ) . 

Contrary to expectations (|Contopoulos et alJll999h we also find that the 
location of the Y-point is inside the Light Cylinder. This is due to the presence 
of a finite pressure plasma at the base of the star. In this case the extent of 
the closed zone is limited by the pressure equilibrium at the Y-point. This 
implies that even for large a the asymptotic wind structure is close to the force- 
free monopole solution, the spindown time can be shorter, and the total energy 
losses larger that what is inferred by applying the dipole formula. 



3. Pulsar Wind Nebulae 

A way to constraint the properties of the wind at large distances is to study 
its interaction with the environment and especially with the surrounding SNR. 
The pulsar wind blows a cavity inside the SNR, filled with a hot relativistic pair 
plasma and magnetic field, known as Pulsar Wind Nebula, which shines in non 
thermal synchrotron emission from radio to X-rays. Recent X-ray images of the 
central region in PWNe have shown the presence of a well defined axisymmetric 
structure known as "jet-torus". A main emission torus is observed in the equa- 
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Figure 3. Right: simulation of the flow structure inside a PWN. The oblate 
shape of the wind region is due to the wind energy distribution. A fast equa- 
torial flow channel survives because magnetization vanishes at the equator. 
At higher latitudes th e flow is collimated tow ard the axis by hoop stresses 
and forms a jet. From lDel Zanna et alJ (2004'). Left: synchrotron map based 
on simulations. Notice the presence of a c entral "knot" an i nner a rc and an 
outer torus. The jet is also visible. From iDel Zanna et alJ |2006). Axes in 
both pictures are in light years. 



torial plane, with the presence of an inner ring of emission; a perpendicular jet 
is present. Sometimes an enhancement of emission close to the pulsar location, 
known as the "knot" is also observed. The observed doppler boosting suggests 
velocities in the torus ~ 0.3 — 0.4c and in the jet ~ 0.6 — 0.7c. Such struc- 
ture is too complex to be interpreted in term of existing ID models. Moreover 
the presence of a jet seems to contradict the result on the wind structure that 
collimation is not efficient. 

Recent numerical simulation (Del Za nna et a has shown that 

such structure naturally arises as a consequence of an anisotropic energy dis- 
tribution and sufficiently high magnetic field in the nebula. Assuming a split 
monopole for the wind energy and lorentz factor distribution, as derived from 
2D simulation of relativistic winds, and a magnetization higher than a ~ 0.01, 
simulated synchrotron maps are able to recover qualitatively all the main fea- 
tures of the jet-torus structure. The torus and inner ring as well as the knot 
are explained in term of region where the flow is moving at a relativistic speed 
toward the observer, while the formation of the jet is associated with collimation 
by hoop stresses inside the nebula. In all simulations a 2D spherical r — 9 grid 
was used, with 100 uniform cells in 0, and a logarithmic scale with 200 cells per 
decade in r. A wind with energy distribution typical of a monopole is injected 
at the base, in an environment representative of SNR ejecta. 

Results suggest that the emission features observed in the nebulae might 
be compared with simulations to constraint the wind properties. Preliminary 
results show that a better agreement is achieved if an unmagnetized region in 
the equatorial plain is assumed. Such unmagnetized region could be due to 
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dissipation of the current sheet (|Kirk fc Skiaeraase which in turn could 

provide a mechanism for accelerating the wind and solving the a problem. 



4. Conclusions 

Recent numerical results have confirmed theoretical predictions about the accel- 
eration of relativistic pulsar winds. The conversion of magnetic to kinetic energy 
appears to be minimal, the flow closely resembles the exact monopole solution, 
there is no evidence for the formation of a collimated energetic jet. The use 
of a more realistic dipolar field does not seem to change significantly the wind 
properties far from the light cylinder. On the contrary the energy and angular 
momentum losses seem to depend only on the integrated value of a, and for 
a > 1 on the amount of open magnetic flux. 

An alternative way to investigate the wind properties at least at large dis- 
tances comes from the study of its interaction with the surrounding SNR. Mod- 
eling of the "jet-torus" structure in PWNe have shown that such features are 
due to the details of the nebular flow, which in turn depends on the wind prop- 
erties. Simulations show that an energy and magnetic field distribution in the 
wind resembling the monopole solution, can indeed explain the main details of 
the observations. 
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